arXiv:1502.07086v3 [quant-ph] 17 Mar 2016 


J. Stat. Mech. (2016) 033102 arXiv: 1502.07086 


Quantum entropy of systems described by non-Hermitian Hamiltonians 

Alessandro Sergi 1,2 >[*] and Konstantin G. Zloshchastiev 3, 0 

1 Dipartimento di Fisica e Scienze della Terra, Universita degli Studi di Messina, Contrada Papardo, 98166 Messina, Italy 

2 School of Chemistry and Physics, University of KwaZulu-Natal, 

Private Bag X01, Scottsville, Pietermaritzburg 3209, South Africa 
3 Institute of Systems Science, Durban University of Technology, P. O. Box 1334, Durban 4000, South Africa 
(Dated: Received: 25 February 2015 [arXiv], 1 September 2015 [JStat]) 

We study the quantum entropy of systems that are described by general non-Hermitian Hamil¬ 
tonians, including those which can model the effects of sinks or sources. We generalize the von 
Neumann entropy to the non- Hermitian case and find that one needs both the normalized and 
non-normalized density operators in order to properly describe irreversible processes. It turns out 
that such a generalization monitors the onset of disorder in quantum dissipative systems. We give 
arguments for why one can consider the generalized entropy as the informational entropy describing 
the flow of information between the system and the bath. We illustrate the theory by explicitly 
studying few simple models, including tunneling systems with two energy levels and non-Hermitian 
detuning. 
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1. INTRODUCTION 

One of the most intriguing problems of statistical me¬ 
chanics is provided by the fact that Hamiltonian re¬ 
versible dynamics is not able to predict any increase of 
the fine-grained entropy, as it would be required by the 
second law of thermodynamics M- However, it has 
been shown that for classical systems such an increase can 
be described through the adoption of non-Hamiltonian 
dynamics 0-[8] with phase space compressibility [§, [nj • 

The difficulties with the reconciliation of the fine¬ 
grained entropy and thermodynamics remain unchanged 
when passing to the realm of quantum mechanics. Here, 
we consider the quantum dynamics originating from gen¬ 
eral non-Hermitian Hamiltonians (NH), known as the 
non-Hermitian approach. This approach is often invoked 
in order to describe quantum systems coupled to sinks 
or sources and it may arise in a variety of contexts, for 
instance, when studying optical waveguides [llj, H3. Fes- 
hbach resonances and par ticles’ disintegration jl3Hl8j | . 
multiphoton ionization 11 Ifl l22l |. and open quantum sys¬ 
tems [23H36I ]. In all such cases, the probability does not 
have to be conserved, in general. 

In this work, we show that the production of the fine¬ 
grained entropy can be naturally predicted within the 
framework of the non-Hermitian approach. In particu¬ 
lar, we extend the definition of the Gibbs-von Neumann 
entropy [371 to the case of systems with non-Hermitian 
Hamiltonians and introduce a “non-Hermitian” entropy 
combining the normalized and non-normalized density 
matrix. In order to illustrate the theory, we explicitly 
consider the analytical solution of some models of interest 
for quantum dynamics. Depending on the model studied, 
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we find that the non-Hermitian entropy can provide the 
expected behavior at large times. 

The structure of this paper is as follows. In Sec. [2] 
we give a brief outline of the density operator approach 
for NH systems. In Sec. [3] we introduce a generalization 
of the Gibbs-von-Neumann entropy that is suitable for 
NH systems, and discuss its features. In Sec. [I] we study 
the NH dynamics of a two-level system in order to illus¬ 
trate the formalism. The discussion of the results and 
the conclusions are presented in Sec. [5] 

2. QUANTUM DYNAMICS WITH 
NON-HERMITIAN HAMILTONIANS 

In the theory of open quantum systems the non- 
Hermitian approach has recently acquired a strong pop¬ 
ularity since it has a different range of applicability from 
the approach based on the Lindblad master equation [38| . 
In order to sketch how the approach unfolds, one can con¬ 
sider a total non-Hermitian Hamiltonian 

H = H - if , (1) 

where both H and F are Hermitian (f is often called the 
decay rate operator), the Schrodinger equations for the 
quantum states |\E f ) and (\H| are written as 


W = 


(2) 

dt{* \ = 


(3) 

Upon introducing a non-normalized density matrix 




(4) 


k 


where Vk are the probabilities of the states (| v I' fc ), (H ,fe |) 
that are compatible with the macroscopic constraints 
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obeyed by the system, the dynamics can be recast in 
terms of the equation 


d t Cl = 




(5) 


where [, ] and {,} denote the commutator and anticom¬ 
mutator, respectively. In the context of theory of open 
quantum systems, the evolution equation for the density 
operator Cl effectively describes the original subsystem 
(with Hamiltonian H) together with the effect of envi¬ 
ronment (represented by T). 

Upon taking the trace of both sides of Eq. © , one 
obtains an evolution equation for the trace of 

<9 t Tr Cl = -jTi (fClj . (6) 


Equation © shows that NH dynamics does not conserve 
the probability. 

As suggested in Ref. (36], one is then led to the intro¬ 
duction of a normalized density matrix, defined as 


P = 


ci 

Tr Cl ’ 


(7) 


that can be used in the calculation of quantum statistical 
averages of arbitrary operators %: ( x)t = Tr ( xp ) • 

As a result of the definition given in Eq. © and the 
evolution equations in © and ©, the normalized density 
matrix obeys a dynamics ruled by the following equation: 


dtp = 



^ {f,p} + ^/5Tr(fp). (8) 


This equation effectively describes the evolution of orig¬ 
inal subsystem (with Hamiltonian H) together with the 
effect of environment (represented by T) and the addi¬ 
tional term that restores the overall probability’s conser¬ 
vation. One can see that due to the last term in this 
equation, the dynamics of the normalized density matrix 
p is nonlinear. A similar nonlinearity was found in the 
evolution equation for the operator averages [39]. More¬ 
over, the appearance of nonlinearities in NH-related the¬ 
ories has also been suggested in Ref. [io|, on the grounds 
of the Feshbach-Fano projection formalism. 

The density operator p, determined by the solution 
of Eq. ©, is bounded and allows one to maintain a 
probabilistic interpretation of the statistical averages of 
operators under non-Hermitian dynamics. Nevertheless, 
the gain or loss of probability associated with the cou¬ 
pling to sinks or sources are properly described by the 
non-normalized density operator Cl. Hence, it turns out 
that one must use both Cl and p in the formalism, one 
to describe gain or loss of probability and the other to 
calculate av erag es. We have already verified in our pre¬ 
vious work [4l| on time correlation functions the need 
to consider both Cl and p in the definitions of statistical 
properties. 


3. QUANTUM ENTROPY 


It is well-known that in the Hermitian case the quan¬ 
tum dynamics is unitary and defined in terms of a nor¬ 
malized density matrix 2 obeying the quantum Liouvillc 
equation of motion: 


d t Z = 



(9) 


The quantum entropy can be defined as 

5 v n =-fcsTr(HlnS), (10) 

where k B is Boltzmann’s constant [37[. The rate of en¬ 
tropy production, derived from the quantum Liouvillc 
equation in ©, is 


d t S vN = -k B Tr (d t z. Ins) = 0. (11) 


While the von Neumann entropy in Eq. m is fit to rep¬ 
resent the properties of equilibrium quantum systems, 
Eq. (fill) implies that the use of the entropy in Eq. (flOl) is 
somewhat more problematic in nonequilibrium dynam¬ 
ics. In fact, in order to agree with the entropy increase 
required by the second law of thermodynamics, one must 
resort to modified definitions of entropy, such as those 
implied by coarse-graining (see Ref. (42J , for example) or 
by the adoption of relevant definitions of entropy |43| . 
This has even led some authors UM to invoke more 
general structures than Hamiltonian ones in order 
to define the microscopic dynamics of statistical systems. 

When a quantum system is coupled to sink or sources, 
NH dynamics can be used. In such a case the straight¬ 
forward adoption of the von Neumann entropy m leads 
to 


S'vn = ~k B (In p) = ~k B Tr (plnp). (12) 

Equation m clearly becomes identical to Eq. m when 
f — > 0 so that W —> H and the dynamics becomes uni¬ 
tary. 

If one takes the time derivative of Eq. ED. uses the 
evolution equation in © and the properties of the trace, 
the following equation for the rate of entropy production 
is obtained: 

d t S vN = —Tr (f p Inp) + |Tr (fp) S vN . (13) 

This equation shows that the non-unitary evolution given 
by Eq. © provides, in general, a non-zero entropy pro¬ 
duction. 

Interestingly, in agreement with our discussion about 
the important role of both Cl and p in non-Hermitian 
dynamics, done in Sec. [2j it is also possible to define the 
entropy as the statistical average of the logarithm of the 
non-normalized density operator: 

Snh = -fcs(lnfi) = -fcsTr(plnQ) = - k B 

Tr Q 

( 14 ) 
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with the evolution of fi naturally given by the linear equa¬ 
tion ©■ One can expect Svn n °t t° be able to catch 
properly the gain or loss of probability because of its sole 
reliance on the bounded p with its nonlinear corrections. 
Instead, the operator In Q can be expected to monitor 
properly the probability evolution. The rate of change of 
S'nh is easily found to be 

d t S NB = —jr—Tr (fpln^ + ^Tr (f>) S NH +2^Tr (fp) 

(15) 

The two entropies are related by the formula 

S'nh = S v n - k B In (Vr , (16) 

therefore, the difference between Snh and S v n is a mea¬ 
sure of deviation of Tr from unity. 

Another important property of the entropy Snh is 
that, unlike the von Neumann entropy m, it is not in¬ 
variant under the complex constant shifts of the Hamil¬ 
tonian that preserve the form of the evolution equation 
for the normalized density operator ©. These constant 
shifts of the Hamiltonian can be regarded as a kind of 
“gauge” transformation, see Appendix. Indeed, if one 
adds to the T operator a constant term that is propor¬ 
tional to the unity operator /, then both the normalized 
density operator 0 and the von Neumann entropy (1121) 
are unchanged; however, the NH entropy acquires a shift 
in terms of a linear function of time: 


r —> r + — hai => 


S V N ~t SvNj 
Snh -t Snh + ks<xt, 


(17) 


where a is an arbitrary real constant, see Appendix for 
details. Such a property can facilitate the computing 
of S'nh for some systems. It also throws light on the 
fact that the NH entropy “remembers” the effects of the 
complex constant shifts, H —» H + cqI (where cq is an 
arbitrary complex number). 


4. EXAMPLES 

In order to demonstrate the behavior of the above- 
mentioned types of entropy, in this section we consider 
few simple models. 


4.1. Models with a constant f operator 

This is the class of models where the Hermitian part 
of the Hamiltonian can be any physically admissible self- 
adjoint operator H whereas the T operator is propor¬ 
tional to the identity operator: 

f = ^70 I, ( 18 ) 

where the parameter 70 is assumed to be real-valued. 


In such models the value of the parameter 70 does not 
affect the time evolution of the normalized density oper¬ 
ator 0. Indeed, Eq. 0 becomes just the conventional 
quantum Liouville equation: 


dtp = 



(19) 


However, equation © reveals that the operator (fT51) does 
affect the evolution of the operator H and of the entropy 
S'nh- Equations © and m become, respectively: 


<9 t TrH = — 7 0 Trf 2 , (20) 

d t S vN = fc B 7 oTr (plnp) + 7 0 S v n = 0. (21) 

Imposing the initial conditions TrS2(0) = 1 and S v n(0) = 
Sy§ = const, and also using the relation (fTCl) . we obtain 


Tr Cl(t) = exp (— 70 1), 

( 22 ) 

S vN (t) = = const, 

(23) 

<Snh(*) = S^n + k B 7ot- 

(24) 


One can see that at large times the trace of Q either di¬ 
verges (at negative values of 70 ) or vanishes (at positive 
values of 70 ) but the conventional von Neumann entropy 
S v n does not reflect this behaviour in any way. On the 
other hand, the entropy S’nh provides more information 
in this regard. For instance, at positive values of 70 the 
trace of H goes asymptotically to zero, describing the 
damping of the probability. In such a case, S'nh grows 
linearly with time as the thermodynamic entropy is ex¬ 
pected to do. Further discussion of these features is given 
in the concluding section. 


4.2. Two-level tunneling model with 
non-Hermitian detuning 

Let us consider a two-level model specified by the Her¬ 
mitian Hamiltonian 


H = —HAa x (25) 

and the T operator 

f = tvy(r z , (26) 

with the total Hamiltonian given by Eq. 0. The pa¬ 
rameters A and 7 are real-valued, with A being also 
positive. The symbols a x and a z denote the Pauli 
matrices: a x = (?q), 07 = (J^). This model is 
the non-Hermitian analogue of the well-known tunneling 
model with detuning phi ] , which finds applications in the 
pseudo-Hermitian and PT-symmetric quantum mechanic 
[45l . l46p . Such kind of models are often used in order 
to effectively describe the dissipative and measurement- 
related phenomena in open quantum-optical and spin 
systems, such as the direct photodetection of a driven 
TLS interacting with the electromagnetic field [H]. 
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where 0 ^ p ^ 1 is a free parameter. It is easy to check 
that this state is pure at p = 0 , 1 and mixed otherwise. 
Solving the evolution equations ©, with the initial con¬ 
dition m, we obtain the following expression for the 
non-normalized density operator 


n = 


2p 2 


2f y {t)a y + f z {t)a z + F(t)I 


(28) 


where we denoted: 


f y [t) = sinh (/it) [ppcosh (pr) — 7 sinh (pr)] , (29) 
f z ( t ) = p [pp cosh (2 pr) — 7 sinh (2 pr)] , (30) 

F(t) = 7 2 cosh (2/rr) — pfry sinh (2 pr) — 1, (31) 


where r = At, 7 = 7 /A, p = 2p — 1 , and the value 
p = 1 / 7 2 — 1 is assumed to be positive throughout the 
paper. One can see that for the chosen initial state m, 
Tr 0 is invariant under the simultaneous transformation 
7 —> —7 and p —> 1 — p, which will manifest itself in the 
behavior of £nh below. Consequently, the normalized 
density matrix is given by: 


- Ml- fz{t) . i f 
P F(t) y+ 2F(t) z+ 2 1 ' 


(32) 


FIG. 1: Entropy S v n (l33l) versus r for the values of 7 = —2 
(top plot) and 7 = 2 (bottom plot), at the different values of 
p: 0.01 (solid line), 1/4 (dashed line), 1/2 (dash-dotted line), 
3/4 (dotted line) and 0.99 (dash-double-dotted line). 


The von Neumann entropy can be computed directly 
from the definition (1121) . It is given by (in units where 
ks = 1 ): 


As an initial state we choose the superposition of the 
ground and excited states 


SvN = -F| 2) (t) In4 1} (*) - F- ] (i) In F^ 1] (t), (33) 


P( 0) = ^(0) = p\e) (e| + (1 - p) | g) (g\ 



where we denoted F±\t) 
1 , 2 , and 


I i±vm/p(t) 


i = 


Fi(t) = 4 sinh 2 (/it) [pp cosh (pr) — 7 sinh (pr )] 2 + p 2 [pp cosh (2/ir) — 7 sinh (2/ir)] 2 , (34) 

F 2 (t) = pp 7 [2 sinh (2 pr) — 7 2 sinh (4pr)] + iy 2 [(p 2 (p 2 + 1) + 1) cosh (4pr) — 4 cosh (2pr)] 

+ \[p i i.P 2 -F)-p 2 {p 2 - 2 )+ 3 ]. ( 35 ) 


The typical profiles of the entropy (l33l) for the initial state 
are shown in Fig. [I] One can see that the entropy 
jS v n tends to zero at large times, regardless of the sign of 

7 - 


The NH entropy, defined in Eq. ED. can be computed 
using the relation (1161) . It turns out to be (in units where 
ks = 1 ): 


<5nh = SVn ~ In \F{t)/p 2 ] 


(36) 
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where S v n is given by (1531) . One can see that if the von 
Neumann entropy remains finite at large times then the 
asymptotical properties of Snh are determined by the 
behavior of F(t), i.e., 

lim Snh = ln/r 2 — lim In F(t) oc — 2/ir, (37) 

£—>•00 £—>oo 

such that Snh tends to a linear function at large times. 




FIG. 2: Entropy Snh (1361) versus r for the values of 7 = —2 
(top plot) and 7 = 2 (bottom plot), at the different values of 
p: 0.01 (solid line), 1/4 (dashed line), 1/2 (dash-dotted line), 
3/4 (dotted line) and 0.99 (dash-double-dotted line). 


The profiles of the entropy (l36l) for the initial state 
(Em are shown in Fig. [2] One can see that the entropy 
decreases with time, regardless of the sign of 7 , as one 
can expect from equations m and (1361) . It is instruc¬ 
tive to compare this with the plots for the trace of the 
non-normalized density matrix, shown in Fig. [3j which 
indicate the flow rate of probability to/from the system. 
For this model, the entropy Snh takes negative values 
at large times. This will be discussed in details in the 
concluding section. 




FIG. 3: Profiles of Trff, where Cl is given by (1281 . versus r 
for the values of 7 = —2 (top plot) and 7 = 2 (bottom plot), 
at the different values of p: 0.01 (solid line), 1/4 (dashed 
line), 1/2 (dash-dotted line), 3/4 (dotted line) and 0.99 (dash- 
double-dotted line). 


4.3. Two-level tunneling model with non-Hermitian 
detuning and asymptotically constant NH entropy 

For the previous two-level model we found out that the 
NH entropy goes to negative values during time evolu¬ 
tion. Here we illustrate that this behavior can be changed 
just by adding a constant decay operator to the F oper¬ 
ator, as explained in the last paragraphs of section [3] 
Thus, the Hermitian part of the model is given by (ESI) 
whereas the F operator is (in units ks = 1): 

F = AycU + h/jl. (38) 

Unless otherwise specified, here and in the following, we 
assume the notation of section ITT21 The initial conditions 
for the evolution equation remain m- 

Using the results given in the Appendix, one can easily 
show that, for this model, both the normalized density 
p and von Neumann entropy S v n are the same as those 
obtained for the model in Sec. roi However, the operator 
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FIG. 4: Entropy Snh (l40l) versus r for the values of 7 = — 2 
(top plot) and 7 = 2 (bottom plot), at the different values of 
p: 0.01 (solid line), 1/4 (dashed line), 1/2 (dash-dotted line), 
3/4 (dotted line) and 0.99 (dash-double-dotted line). 


FIG. 5: Profiles of Tr ft, where ft is given by HO) , versus t, 
for the values of 7 = — 2 (top plot) and 7 = 2 (bottom plot), 
at the different values of p\ 0.01 (solid line), 1/4 (dashed 
line), 1/2 (dash-dotted line), 3/4 (dotted line) and 0.99 (dash- 
double-dotted line). 


(l28l) acquires the factor exp (— 2 /it): 


3-2 Vt 


2/i 2 


2 fy(t)&y + fz(t)a z + F(t)I 


(39) 


Therefore, using the transformations in Eq. m, we ob¬ 
tain 


-Snh = SVn - In [F(f)//i 2 ] + 2/ir, (40) 

where .SVn is given by (1531) . The asymptotical value is 
given by 


lim S'nh = In 

t—y 00 


2/i 2 

7(7 - fip)_ 


(41) 


Unlike its analogue in Sec. 14.21 the S'nh entropy for this 
model is bounded. As a matter of fact, considering the 
more general form of the T operator, T = /ry<r 2 + khfil 
(where k is an arbitrary real number), one can establish 
that the value k = 1 acts as a critical threshold: the NH 
entropy decreases asymptotically if k < 1 (see the model 
in Sec. l4~2ll . while it increases if k > 1 (see the model in 
Sec. 14.41) . 


The profiles of the von Neumann entropy, the NH en¬ 
tropy and the trace of the operator fl for this model are 
shown in Figs. [T][I]and[5l respectively. On can see that, 
as for the model in Sec. 14.21 the entropy S'nh may takes 
negative values for a certain range of parameters (when 
the trace of Cl goes above one). This will be discussed in 
details in the concluding section. 


4.4. Two-level tunneling model with 
non-Hermitian detuning and increasing NH entropy 

In the two previous cases, we found that the NH en¬ 
tropy can assume a negative infinite value or a constant 
value. Here we consider another model that, instead, pro¬ 
vides an entropy that asymptotically increases with time. 
Also in this case, such a behavior can achieved just by 
adding a constant term to the T operator, as explained 
in the last paragraphs of section [3] 

The Hermitian part of this model is still given by Eq. 
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FIG. 6 : Entropy Snh (l44l) versus r for the values of 7 = — 2 
(top plot) and 7 = 2 (bottom plot), at the different values of 
p: 0.01 (solid line), 1/4 (dashed line), 1/2 (dash-dotted line), 
3/4 (dotted line) and 0.99 (dash-double-dotted line). 


(I?5l) whereas the f operator is (in units ks = 1): 


r = 


h'yd'z + —hpi. 


(42) 


As one can see, the constant term multiplying the iden¬ 
tity operator in the definition of F above is | hp. Al¬ 
though, this case is qualitatively similar to any other one 
with the constant term’s coefficient larger than hp. The 
initial conditions for the evolution equation remain (1271) . 

Using the results given in the Appendix, one can easily 
show that for this model both the normalized density p 
and von Neumann entropy S v n are the same as for the 
models in Secs. 14.21 and roi whereas the operator (1281) 
acquires the factor exp (—3 pt): 


2p 2 


.,-3/it 


2f y (t)a y + f z {t)a z + F(t)I 


Therefore, using El. we obtain 

Snh = S'vN ^ In [F(t)/p 2 ] + 3/ir, 


(43) 


(44) 


where S'vN is given by (l33l) . The asymptotical value is 
given by 

lim Snh = ln/r 2 — lim In F(t) + 3 /it oc pr , (45) 

£—>■+00 £—>• 00 

such that Snh tends to a linear function with a positive 
coefficient at large times. 

The profiles of the von Neumann entropy, NH entropy 
and the trace of the operator 12 for this model are shown 
in Figs. CO© and[3 respectively. 




FIG. 7: Profiles of Trl2, where Cl is given by 03 , versus t, 
for the values of 7 = —2 (top plot) and 7 = 2 (bottom plot), 
at the different values of p: 0.01 (solid line), 1/4 (dashed 
line), 1/2 (dash-dotted line), 3/4 (dotted line) and 0.99 (dash- 
double-dotted line). 


5. DISCUSSION AND CONCLUSIONS 

In this paper we have provided a generalized formu¬ 
lation of the quantum fine-grained entropy for systems 
described by non-Hermitian Hamiltonians. We have 
adopted a straightforward generalization of the von Neu¬ 
mann entropy, defined in terms of the normalized den¬ 
sity matrix (obeying a nonlinear equation of motion), 
and introduced another definition of an entropy, Snh, 
in terms of the normalized average of the logarithm of 





























the non-normalized density matrix. We have shown that 
in both cases the entropy production is non-zero. How¬ 
ever, we have found that it is Snh that properly captures 
the physical behavior of the probability and disorder in a 
system in presence of sinks or sources described by non- 
Hermitian Hamiltonians. 

In Sec. [4] we have studied some models in order to 
illustrate the different behavior of the S v n and Snh en¬ 
tropies. In particular, for the models in Secs. 14.2114.41 
both the normalized density operator, p, and the von 
Neumann entropy, S v n, do not change while the entropy 
Snh does. We show that at large times the value Snh 
decreases asymptotically for the model of Sec. 14.21 it 
tends to a finite value for the model in Sec. m and it 
increases for the model in Sec. 14.41 

The results of the present work, when considered to¬ 
gether with our previous studies [33, Hg, EH , allow us to 
draw a certain number of conclusions. Non-Hermitian 
dynamics is able to describe in the quantum realm the 
production of entropy in a way similar to what non- 
Hamiltonian dynamics with phase space compressibility 
does in the classical realm. Non-Hermitian dynamics 
also seems to need both the normalized density opera¬ 
tor, p, and non-normalized one, Cl, in order to provide 
a proper statistical theory. The non-normalized density 
operator, defined as a solution of Eq. ©, captures some 
important features of the decay process, such as the non¬ 
conservation of probability in the (sub)system and its 
“leakage” into the surrounding environment. The nor¬ 
malized density operator guarantees that the probabilis¬ 
tic interpretation of averages can be maintained. In this 
regard, the entropy Snh combines both operators in a 
proper way and can signal the expected thermodynamic 
behavior of an open system. The entropy Snh also seems 
more suitable for describing the gain-loss processes that 
are related to the probability’s non-conservation, since 
it contains information not only about the conventional 
von Neumann entropy S v n but also about the trace of 
the operator Cl, according to the relation (1161) . Assum¬ 
ing that S v n is bound at large times, the NH entropy 
grows when Tr Q decreases, also it takes positive values 
if Tr Cl < 1 and negative ones otherwise. Hence, one can 
say that Snh describes the flow of information between 
the system and the bath. 

Further studies are needed to understand whether and 
when Snh may deserve a complete quantitative thermo¬ 
dynamic status as well as whether there might be viable 
physical interpretations of a negative entropy (not nec¬ 
essarily given in terms of the number of occupied micro¬ 
scopic states). 
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Appendix: Hamiltonian shift transformations and 
entropy 

Following the discussions presented in Refs. [351, H], 
let us consider the following transformation of the T op¬ 
erator 

f = f' + i Hal, (A.l) 

where a is an arbitrary real constant and I is the unity 
operator. This transformation is a subset of the trans¬ 
formation 

H = U' + c 0 I, (A.2) 

Co being an arbitrary complex number, which is the non- 
Hermitian generalization of the energy shift in conven¬ 
tional quantum mechanics. Therefore, in Refs. [35j, [3fJ 
it was called the “gauge” transformation of the Hamilto¬ 
nian 0, whereas the terms of the type cqI can be called 
the “gauge” terms. 

In Ref. [3(J it was shown that the equation (0 is in¬ 
variant under the transformation (IA.1I) . therefore, one 
immediately obtains 

P = p ', Svn = S(, N , (A.3) 

therefore the von Neumann entropy is not affected by 
the transformation (IA.1I) . One can see that any infor¬ 
mation regarding the shifting of the total non-Hermitian 
Hamiltonian is lost if one deals solely with the normalized 
density operator. 

However, the evolution equation ([6]) is not invariant 
under the shift (IA.1I) . If T~L is time-independent then, 
substituting JA3D into ©> we obtain that the non- 
normalized density acquires an exponential factor: 

Cl = Cl' exp (—at), (A.4) 

therefore, recalling the relation (nsi, we obtain 

Snh = SvN ~ ks In Tr Cl = Sn H + ksat, (A.5) 

which indicates that any lost information about the shift¬ 
ing term in Eq. (IA.2I) in the total non-Hermitian Hamil¬ 
tonian, due the normalization procedure in Eq. 0 , can 
be recovered by means of the NH entropy. 
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